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Introduction 


This  program  grew  out  of  previous  work  aimed  at  understanding  and  producing 
materials  with  large  optical  constants  (dispersion  and  loss)  in  the  electromagnetic  micro- 
wave  region.  Ferro  electrics  showed  promise  in  this  area,  particularly  if  their  paraelectric- 
ferroelectric  transition  temperatures  lie  in  the  ambient  region.  We  were  able  to  explain  this 
behavior  as  arising  from  ionic  motions  in  the  double-well  potential  which  we  now  believe  to 
be  intrinsic  to  virtually  all  structural  transformations  in  ionic  insulators.  However,  only  if 
this  motion  has  a  uniform  polar  component,  i.e.,  is  incipiently  ferroelectric,  can  it  couple 
directly  to  external  electromagnetic  radiation. 

The  basic  thrust  of  the  new  program  was  to  attempt  to  predict  and  explain  ferroelectric 
behavior  from  first  principles  based  primarily  on  the  use  of  the  local  density  approximation 
of  which  we  employed  the  simplest,  the  Gordon-Kim  modified  electron  gas  technique. 

Novel  Halide-Based  Systems 

While  the  primary  interest  was  in  known  ferroelectrics,  i.e.,  oxide-based  perovskites 
such  as  BaTiOa,  our  intial  focus  was  on  novel  halide-based  pervoskites,  one  of  which, 
NaCaFs,  we  had  identified  theoretically  by  employing  Gordon-Kim  potentials  and  molecular 
dynamics,  to  be  a  strong  primary  ferroelectric  below  ~  600K.  Since  the  pure  compound  has 
yet  to  be  synthesized  by  standard  techniques,  we  first  began  to  study  binary  mixed  fluorides 
in  which  NaCaFs  is  one  component  and  the  other  is  a  fluoride  which  can  be  synthesized  in 
its  pnre  form:  specifically,  we  chose  first  KCaFs  and  then  NaMgFs. 


2 


The  resultant  studies  are  described  in  the  following  pair  of  abstracts: 
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THE  EFFECT  OF  K  DEFECT  CLUSTERS  ON  THE 
FERROELECTRIC  PHASE  TRANSITIONS  IN  NaCaF3 

Z.  MO,  J.  W.  FLOCKEN,  W.  N.  MEI  and  R.  A.  GUENTHER 

Unh-ersiT)^  of  Nebraska,  Omaha,  Nebraska,  661 S2  U.S.A. 

and 

J.  R.  HARDY 

Universin  of  Nebraska,  Lincoln,  Nebraska,  6S5S8-0I11,  U.S.A. 

and 

L.  L.  BOYER 

Naval  Research  Laboraioiy,  Washington  D.C.,  20357-5000  U.S.A. 

(Received  March  6,  1991;  in  final  form  June  29,  1991) 

We  have  performed  a  number  of  molecular  dynamics  computations  on  -20  ion  samjples  of  NaCaFj 
containing  K  substitutional  defects  in  various  concentrations  and  conSguraiions.  Ab  iniiio  potentials 
with  DO  adjustable  paramieters  were  used  to  obtain  the  short  range  interactions  between  ion  pairs.  Our 
results  indicate  that  the  presence  of  these  defects,  in  most  cases  has  a  profound  effect  on  the  charac¬ 
teristics  of  the  polar  state,  including  a  diminished  sharpness  of  the  transition  temperature.  We  believe 
that  these  results  m.ay  lend  so.me  insight  into  the  mechanism.s  which  produce  •'relaxor"  ferroelectrics. 

Keywords:  ferroekcirks,  perovskiies.  molecular  dynamics,  defects 
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DIFFUSIVE  PHASE  TRANSITIONS  IN  MIXED 
HALIDE  FERROELECTRICS 


J.  W.  FLOCKEN  and  Z.  MO 

University  of  Nebraska,  Omaha,  Nebraska  661 82,  USA 

and 

J.  R.  HARDY  and  H.  M.  LU 

University  of  Nebraska,  Lincoln,  NE  68588-0111,  USA 

(Received  March  6.  1992;  in  final  form  June  22,  1992) 

Molecular  dynamics  simtulations  have  been  performed  for  the  Nal'fgi.^CajF;.  m.aierials  where  Ca  ions 
were  substituted  for  Mg  in  NaMgF,  in  vjrying  concentrations  of  randomly  distributed  defects  repre- 
sentirtg  constitutional  crossover.  Sm.all  clusters  of  Ca  defects  in  NaMcF,  and  clusters  of  .Ms  defects  in 
.N'aCaFj  were  also  studied.  First-principles  calculations  were  used  to  calculate  she  potentials  between 
ion  pairs  in  these  samples.  The  change  of  defect  concentration  results  in  a  gradual  shift  from  nonfer- 
roelectric  behavior  of  NaMgF,  to  ferroelectric  behavior  of  NaCaF;..  Doping  in  these  m.aterials  leads  to 
a  "sm, eared  out"  transition  "range"  instead  of  a  sharp  transition  temperature — the  so-called  diffusive 
transition. 

Keywords:  diffusive  phase  iranshion,  perovskhe 
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While  ferroelectric  halide  perovskites  have  yet  to  be  fabricated,  other  systems,  notably 
BaMgF4,  do  show  polar  phases;  moreover,  the  polarity  is  “switchable”  by  an  external  field. 

We,  therefore,  initiated  and  carried  through  a  fairly  extensive  study  on  these  systems 
which  is  described  in  the  following  pair  of  abstracts: 
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MOLECULAR  DYN.AMICS  SIMULATION  OF  BaMgF4  and  Ba]vlnF4 


Z.  MO  AND  J.  W.  FLOCKEN 

University  of  Nebraska,  Omaha,  Nebraska,  68182  USA 
DORIAN  M.  H.ATCH,  Brigham  Young  University,  Provo,  Utah, 
84602  USA 
J.  R.  HARDY 

University  of  Nebraska,  Lincoln,  Nebraska,  68588-0111,  USA 

(Received  August  9,  199S) 

Abstract  We  have  performed  molecular  dynamics  simulations 
on  BaMgF4  and  BaMnF4,  both  of  which  have  been  reported  to 
have  polar  structures,  with  Tc  above  the  melting  point.  We 
have  determined  the  ground  state  structure  and  verified  the 
polarization  mechanism  described  by  Keve  e:  alA.  Polarizations 
of  8.5  uC/cm2  are  found  for  the  Mg  compound  at 
6  K  and  11.7  jiC/cm^  for  the  Mn  compound.  We  find  no 
indication  of  ferroelectric  phase  transitions  between  0  K  and 
1200  K  for  either  material.  We  are  able  to  reverse  the 
polarization  of  the  Mg  compound  by  simulating  a  external 
field,  while  the  Mn  based  material  distorts,  but  does  not 
reverse  polarity.  We  have  simulated  substitutions  of  Be  for  Mg 
and  Mn  in  these  materials.  We  report  on  these  results. 
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Molecular-dynamics  simulations  of  some  BaJiT4  compounds 

John  Flocken,  Z.  Mo,  and  W.  N.  Mei 
Universiry  of  Sebrcska  at  Omaha,  Omaha,  Sebraska  6S182 

J.  R.  Hardy 

Univenity  of  Nebraska  at  Lincoln,  Liricoln,  Nebraska  68588-011] 

Dorian  M.  Hatch 

Brigham  Young  University,  Provo,  Utah  84602 
(Received  21  September  1993) 

We  have  carried  out  molecular-dynamics  simulations  on  compounds,  where  X  is  Mg,  Mn,  or 

Zn.  Ab  initio  potentials,  with  no  adjustable  parameters,  were  used  to  obtain  short-range  interactions  be¬ 
tween  ion  pairs.  We  found  a  polar  ground-state  structure  which  is  in  agreement  with  the  Al^am  space 
group  reported  experimentally.  We  were  able  to  reverse  polarization  in  BaMgF*  at  high  temperatures, 
using  large  fields,  but  were  unable  to  reverse  polarization  in  the  other  compounds.  The  second-order 
phase  transition  in  the  Mn  compound  at  250  K  was  reproduced.  We  believe  this  to  be  the  first  extension 
of  molecular  dymamics  to  materials  consisting  of  chains  of  F  octahedra. 

While  we  have  subsequently,  using  State  of  Nebraska  Research  Initiative  and  funding 
under  the  NSF  EPSCoR  program,  attempted  to  synthesize  some  of  these  mixtures  exper¬ 
imentally,  we  have  found  that  if  Na"*"  is  present,  it  is  segregated  as  NaF  together  with 
similarly  segregated  CaF2. 


Oxide  Ferroelectrics 

In  parallel  with  this  work  on  halides,  we  finalized  and  made  specific  our  work  on  the 
dielectric  properties  of  linear  double  well  systems  by  applying  it  to  SbSI.  This  is  summarized 
in  the  following  abstract: 
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A  DOUBLE  WELL  OSCILLATOR  MODEL  FOR  THE 
FERROELECTRIC  PHASE  TRANSITION  IN  SbSI 

J.  W.  FLOCKEN  and  R.  A.  GUENTHER 
Physics  Department,  University  of  Nebraska  at  Omaha,  Omaha,  Nebraska  6S1S2 

and 

J.  R.  HARDY 

Department  of  Physics,  University  of  Nebraska  at  Lincoln, 

Lincoln,  Nebraska  685S8 

and 

L.  L.  BOYER 

Naval  Research  Laboratory,  Washington,  D.C.  20375 
(Received  December  17,  1991] 

The  dielectric  response  spectrum  for  antimony  sulfoiodide  has  been  computed  using  a  model  in  which 
the  material  is  simulated  by  a  damped  double  well  oscillator.  It  is  found  that  with  the  use  of  sirntple 
dam.ping  and  relaxation  terms  and  realistic  well  parameters  the  experimental  dielectric  response,  in¬ 
cluding  the  central  peak,  can  be  reproduced  quite  closely.  The  choice  of  well  depth,  in  particular,  plays 
a  major  role  in  determining  the  frequency  of  the  first  peak  above  the  soft  mode. 

As  the  work  by  other  groups  funded  under  this  program  progressed,  it  produced  a 
disturbing  anomaly.  While  their  local  density  approximation  LAPW  total  energy  calcula¬ 
tions  could  produce  double- well  ferroelectric  structures  for  BaTiOs  etc.,  this  was  only  true 
if  the  experimental  lattice  constants  were  imposed.  At  the  exact  theoretical  volume,  that 
of  minimum  static  energy,  the  minima  were  either  absent  or  implausibly  shallow. 

This  raised  the  question  of  how  valid  any  calculation  which  imposes  periodic  boundary 
conditions  was  for  a  ferroelectric,  owing  to  the  infinite-range  effects  associated  with  plane 
polar  surfaces. 

We  thus  made  a  specific  study  of  this  question  using  a  system  of  NaNOa  molecular 
ions.  The  results,  described  in  the  abstract  below,  indicate  that  this  is  a  genuine  problem 
(whether  or  not  it  influences  the  double  wells  in  BaTiOa  and  other  oxide  perovskites),  and 
its  understanding  is  of  serious  practical  importance  in  the  context  of  thin-film  ferroelectric 
memories. 
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FERROELECTRIC  TRANSITIONS,  MOLECULAR 
DYNAMICS  AND  INTRINSIC  SURFACE  PROBLEMS 

J.  R.  HARDY  and  H,  M.  LU 

Department  of  Physics  and  Center  for  Electro-Optics,  University  of  Nebraska, 
Lincoln,  Nebraska  68588-0111,  USA 

(Received  March  6,  1992;  in  final  form  June  25,  1992) 

Conventional  supercell  molecular  dynamics,  if  applied  to  a  system  where  one  has  reasonable  potentials 
and  which  is  known  experimentally  to  have  a  ferroelectric  ground  state,  can  fail  to  reproduce  that  state. 
We  shall  show  that  such  behavior  stems  from  a  failure,  intrinsic  in  this  type  of  molecular  dynamics,  to 
handle  the  surface  effects  intrinsic  to  Coulomb  interactions.  We  will  discuss  how  this  limitation  can  be 
addressed  and  its  possible  implication  for  ground  state  energies  calculated  by  using  periodic  boundary 
conditions. 

Keywords:  ferroelectrics,  supercell  MD 


While  we  were  carrying  out  our  studies  on  non-oxide-based  systems,  we  were  continu¬ 
ally  developing,  in  conjunction  with  Dr.  Boyer  at  the  Naval  Research  Laboratory,  extensions 
and  applications  of  the  PIB  (potential  induced  breathing)  model,  which  represented  an  early 
attempt  to  study  oxides  allowing  for  the  major  effects  of  charge  distortions  on  the  po¬ 
tential  energy  surfaces  for  all  oxide-based  systems:  notably,  ferroelectric  perovskites. 

Unfortunately,  this  had  only  been  carried  to  the  point  of  determining  the  effect  of  PIB 
on  rotational  (nonpolar)  instabilities,  and  developing  a  PIB-inclusive  molecular  dynamics 
program  when  funding  was  terminated. 

Three  specific  studies  are  described  in  the  following  abstracts,  the  last  representing 
our  definitive  results  within  the  PIB  limitations: 
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MOLECULAR  DmAMICS  SDvRJLATION  OF  THE  PH.4SE 
TRANSITION  IN  LaAlOs 


J.  W.  FLOCKEN  AND  Z.  MO 

University  of  Nebraska,  Omaha,  Nebraska,  68182  USA 
J.  R.  HARDY 

University  of  Nebraska,  Lincoln,  Nebraska,  68588-0111,  USA 

(Received  August  9,  1993) 

Abstract  We  have  used  molecular  dynamics  and  lattice  statics 
to  study  LaAlOs  in  the  temperature  range  0  -  1500  K  .  Short 
range  pair  potentials  were  obtained  using  the  Gordon-Kim^ 
approach.  The  O'-  ion  charge  density  was  stabilized  using 
Watson  Sphere  potentials  ranging  from  0.7  to  0.8  Har./e. 
within  this  range,  a  phase  transition  is  found  to  occur  at  the 
temperature  reported  experimentally  (800K),  however,  we 
find  the  ground  state  to  have  a  tetragonal  (Pmnb)  structure 
ratherjhan  the  observed  R3C  configuration.  It  is  found  that 
the  R3C  structure  can  be  stabilized  by  transferring  charge 
from  the  0'^  ion  to  the  La'*‘3  ion  in  the  long  range  potential 
suras. 


Ferroelectrics,  1994,  N’ol.  153,  pp.  79-84  ©  1994  Gordon  end  Breach  Science  Publishers  S..A.. 

Reprints  available  directly  from  the  publisher  Printed  in  the  United  States  of  .America 

Photocopying  permitted  by  license  only 


A  STLTDY  OF  THE  STRUCTURE  AND  SUABILITY  OF  SiCu02 


J.W.  FLOCKEN 

University  of  Nebraska,  Omaha,  Nebraska,  68182  US.\ 

J.R.  HARDY 

University  of  Nebraska,  Lincoln,  Nebraska,  68588-0111,  USA 

(Received  August  9,  1995) 

Abstract  We  have  investigated  the  structure  and  stability  of 
SrCu02  using  both  lattice  statics  and  dynamics  techniques  as 
well  as  molecular  dynamics.  Short  range  interatomic  pair 
potentials  were  obtained  using  Gordon-Kim^  techniques.  Lattice 
dynamics  indicates  that  the  infinite  layer,  high  Tc  parent 
structure  of  Azuma  et.  al.  -  should  not  be  stable  for  this 
material,  with  the  largest  instabilities  appearing  in  the  plane 
perpendicular  to  the  Cu-Cu  bonds  (i.e.  in  the  ab  plane).  The 
ground  state  structure  and  temperature  variation  of  physical 
parameters  obtained  from  molecular  dynamics  simulations  will 
be  discussed,  as  will  the  effects  of  varying  ionicity. 
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INTERPLAY  OF  IONIC  BREATHING  AND 
ROTATIONAL  INSTABILITIES  IN  PEROVSKITES* 

H.  Z.  CAO  and  J.  R.  HARDY 

Department  of  Physics,  University  of  h’ebraska-Lincoln, 

Lincoln,  NE  68588-0111,  USA 

(Received  April  I.  1994;  in  final  form  June  20,  1994) 

In  a  previous  work'  we  have  used  parameter-free  poientials.  developed  by  the  Gordon-Kim  technique, 
to  study  rotational  instabilities  in  halides-based  perovskiies  by  molecular  dynamics.  These  calculations 
assumed  rigid  spherical  ions  unchanged  by  incorporation  in  the  crystals.  We  have  now  extended  this 
approach  to  include  the  potential-induced  breathing  (PIB)-  associated  with  purely  radial  distortions  in 
the  ionic  charge  densities  induced  by  the  electrostatic  Madelung  field  in  the  lattice.  In  addition  we  have 
examined  the  effects  of  using  different  forms  of  density  functionals’  in  the  Gordon-Kim  potentials. 
Results  will  be  presented  for  both  halides-based  and  oxygen-based  perovskiies. 


Keywords:  Molecular  dynamics,  Gordon-Kim  potentials,  perovskiies 

However,  Dr.  Boyer  has  since  made  major  progress  with  generalizations  of  the  PIB 
approach,  and  we  are  currently  continuing  this  collaboration  with  local  funding. 
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